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Abstract 

The  acoustical  performance  of  a  submerged  linear  array  of  spherical  shells  in  shallow 
water  is  examined  by  combining  the  T-Matrix  method  of  solving  for  multiple  acoustic 
interactions  among  separate  bodies,  with  a  model  for  the  ocean  as  a  fluid  layer  over  a 
half-space  of  a  distinct  fluid.  The  system  of  source  and  reflected  waves  is  analyzed  by 
standard  approximate  contour  integration  techniques,  valid  for  the  far  field  of  the  array. 
Calculations  using  the  results  show  the  effects  of  shallow  water  on  array  performance. 


INTRODUCTION 

The  prediction  of  the  performance  of  an  array  of  active  transducers  in  a  variety  of  conditions 
is  a  significant  practical  problem,  and  many  methods  have  been  applied  to  this  and  related 
problems. 

Here  we  model  active  array  performance  in  shallow  water  by  coupling  two  representations 
of  the  acoustic  field:  a  T-matrix  representation  near  the  array  and  a  representation  of  the 
water  “sound  channel”  with  a  free  surface  and  a  bottom  of  uniform  depth,  which  overlies  a 
fluid  half-space  of  differing  material  properties,  i.e,  the  Pekeris  model. 

The  applicability  of  the  Pekeris  model  has  been  discussed  at  length  in  the  memoir  of  Prop¬ 
agation  of  Sound  in  the  Ocean  by  Ewing,  Worzel,  and  Pekeris  [1]  and  has  been  instrumental 
in  explaining  acoustical  features  of  sound  propagation  in  coastal  waters.  More  recently,  Tol¬ 
stoy  and  Clay  [2]  have  discussed  more  accurate  models  for  a  shallow  water  channel  which 
include  the  effects  of  attenuation  in  the  shallow  water  “waveguide”  by:  (a)  the  introduction 
of  dissipation  in  the  fluid  half-space  model  for  the  ocean  bottom;  (b)  introduction  of  a  third 
or  multiple  fluid  layers;  or  (c)  allowance  for  the  transmission  of  shear  waves  into  the  bottom 
half-space.  Extensions  of  the  present  work  could  be  applied  to  any  of  these  more  elaborate 
models  with  a  subsequent  increase  in  the  complexity  of  the  problem.  For  the  purpose  of 
demonstrating  the  effects  of  coupling  between  radiating  elements  in  a  shallow  water  array, 
however,  the  simplest  model  (the  Pekeris  model)  has  been  chosen. 
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There  have  been  a  series  of  papers  concerning  sound  propagation  in  shallow  water  pro¬ 
duced  by  an  array  of  projectors  [3, 4, 5, 6]  but  in  each  case,  coupling  between  the  individual 
elements  in  the  array  are  neglected.  In  the  present  work,  this  coupling  is  explicitly  taken 
into  account  by  a  T-matrix  type  analysis. 

The  T-matrix  approach  accounts  explicitly  for  the  mutual  interactions  (including  multiple- 
scattering  effects)  of  elements  from  the  array,  as  well  as  the  individual  radiation  of  each 
element  caused  by  an  internal  driving  force.  A  radiated  acoustic  field  in  terms  of  outgoing 
spherical  harmonic  waves  from  each  element  is  found  which  is  consistent  with  an  array  sub¬ 
merged  in  a  fluid  of  infinite  extent.  The  methodology  is  similar  to  that  reported  by  Scandrett 
and  Canright  [7]  in  which  source  level  calculations  for  arrays  of  spherical  shells  in  an  infinite 
fluid  medium  were  found. 

The  method  of  images  yields  a  superposition  of  reflected  plane  waves  for  each  spherical 
harmonic  excited  by  an  element  of  the  array,  with  the  reflection  coefficients  chosen  to  give  the 
correct  behavior  at  the  free  surface  and  bottom  surface.  Intermediate  steps  taken  to  account 
for  the  shallow  channel  conditions  consist  of:  representing  a  radiated  spherical  harmonic  in 
terms  of  a  superposition  of  plane  waves  over  real  and  complex  angles;  constructing  an  infinite 
set  of  images  in  order  to  satisfy  the  free  surface  and  bottom  surface  boundary  conditions; 
formally  summing  the  series;  and  evaluating  the  acoustic  far  field  by  an  application  of  contour 
integration.  This  sequence  of  steps  has  been  done  in  Brekhovskikh  [8],  and  (using  a  slightly 
different  formulation)  by  Ewing,  Jardetzky,  and  Press  [9],  and  Pekeris  [1],  for  the  case  of  a 
monopole  point  source  in  a  Pekeris  model  for  a  shallow  water  channel. 

It  should  be  pointed  out  that  the  coupled  approach  taken  in  the  current  work  neglects 
the  effects  of  surface  or  bottom  reflections  on  the  spherical  shells  themselves.  Such  effects 
will  be  small  if  the  array  is  located  several  wavelengths  from  the  surface  and  bottom,  but  if 
such  effects  were  deemed  important  (perhaps  for  the  case  of  an  array  near  the  free  surface  of 
the  shallow  water  channel),  the  T-matrix  approach  could  easily  be  extended  to  include  the 
“free  surface  image”  of  near  surface  array  elements. 

FORMULATION 

The  Pekeris  model  for  a  shallow  water  waveguide  is  employed  in  which  the  upper  fluid  layer 
(p  =  1000  kg /m3,  c  =  1524  m/s)  lies  over  the  bottom  fluid  halfspace  (pi  =  2000A;g/m3, 
ci  =  1600  m/s).  The  bottom  is  assumed  to  be  at  a  constant  depth  of  H  =  10 m.  Embedded 
in  the  fluid  layer  is  an  array  of  steel  spherical  thin  shells  which  are  forced  to  oscillate  in  a 
breathing  mode  at  a  frequency  485  Hz  (cc).  The  thin  steel  shells  have  a  radius  of  .5  m(a), 
are  .0026162  m  (h)  thick  and  have  the  material  properties  :  E  =  2.07  x  10 11  N/m2,  v  =  .3, 
and  ps  =  7669  kg/m3. 

The  governing  equations  used  are: 

Fluid: 


V2p  +  k2p  =  0  ,  k  —  —  ,  0  >  z  >  —H 

c 

V2pi  +  k\pi  =  0  ,  h  =  — 

Cl 


2 


z  <  —H 


(1) 

(2) 


(3) 

(4) 

(5) 

(6) 
(7) 


Fluid/Structure  Interface: 

Vp  •  Hj  =  iujpwj  (8) 

where 


and  Hj  is  the  unit  normal  to  the  jth  shell. 

A  periodic  e~lut  time  dependence  is  assumed  in  the  solution,  and  all  distances  are  scaled 
by  the  value  of  the  wavenumber  k,  which  equals  2  for  the  values  of  uj  and  c  given  above, 
resulting  in  ka  =  1  for  the  prescribed  problem,  with  a  scaled  layer  depth  of  twenty.  (See  fig. 
1  for  a  schematic  of  the  problem  as  stated.) 

METHODS 

In  a  previous  work  of  the  authors  [7],  source  level  calculations  of  a  low-frequency  active  array 
of  elastic  spherical  shells  with  specified  internal  forcing  were  found.  In  that  work  however, 
results  were  applicable  only  to  arrays  of  projectors  in  an  infinite  fluid  medium.  The  acoustic 
pressure  held  produced  by  the  array  was  represented  by  a  superposition  of  outgoing  spherical 
waves  for  which  far-held  source  level  calculations  were  made. 


3 


The  starting  point  of  the  present  analysis  begins  with  this  series  representation  of  the 
acoustic  field  for  an  unbounded  fluid  medium. 

L  N  n 

p(x)  =  ee  y  sLyi>(fcn)p,r(coS9I)e”"»‘  (io> 

1=1  n= 0  m=—n 

where  the  harmonic  time  dependence  e~luJt  is  omitted.  The  l  index  refers  to  individual 
elements  of  the  array  (with  the  spherical  coordinates  ( ri,0i,(f>i )  relative  to  the  position  of 
the  Ith  element)  while  the  m,  n  indices  are  over  the  spherical  harmonics  excited  by  the  Ith 
spherical  shell,  and  whose  number  is  truncated  by  the  value  of  N.  The  amplitudes  ( Blnm ) 
incorporate  the  forced  radiation  of  the  Ith  spherical  shell,  as  well  as  reflection  of  acoustic 
pressures  from  the  remaining  spherical  shells  of  the  array. 

The  shallow  water  effects  are  included  by  considering  a  single  spherical  harmonic  excited 
by  a  single  spherical  shell,  with  a  view  to  determine  how  the  acoustic  field  produced  by  this 
single  harmonic  must  be  modified  in  order  to  satisfy  the  free  surface  and  bottom  boundary 
conditions  present  in  the  shallow  water  channel.  Summation  over  all  harmonics  and  all 
spherical  shells  can  then  be  done  after  the  effect  of  a  single  harmonic  of  a  single  array 
element  has  been  found.  The  effect  of  the  mutual  array  element  interactions  is  present  by 
virtue  of  the  manner  in  which  the  Blnm  have  been  calculated. 

While  the  radiated  field  is  expressed  in  terms  of  spherically  spreading  waves,  the  spherical 
waves  must  be  reformulated  in  terms  of  plane  waves  in  order  to  account  for  the  angle  depen¬ 
dence  of  the  reflection  coefficient  at  the  interface  between  two  fluid  media.  Each  outgoing 
spherical  wave  can  be  written  as  a  superposition  of  plane  waves  as  follows  [10,11]: 

j  n  p2n  r  — ioo 

h£\kr)P™ (cos  0)eim*  =  (±l)m+n  —  /  /  eikrcosiP™(cosa)eimP  sina  da  d(3  (11) 

2ti  Jo  Jo 

where  the  double  integral  is  over  plane  waves  in  all  directions  (including  complex  angles, 
corresponding  to  “inhomogeneous”  plane  waves)  and  cos  7  =  sin  d  sin  a  cos(</>— /3)±cos0  cosa 
with  the  upper  signs  for  plane  waves  propagating  above  the  source  (z  >  0),  and  the  lower 
sign  for  plane  waves  propagating  below  the  source  {z  <  0).  (The  contour  integral  in  a  runs 
along  the  real  axis  to  ir/2,  then  parallel  to  the  imaginary  axis.) 

The  reflection  {V\ )  and  transmission  (T)  coefficients  for  a  plane  wave  of  unit  amplitude 
striking  a  planar  interface  between  two  fluids  of  differing  material  properties  at  an  angle  a 
with  respect  to  the  normal  of  the  interface  are  given  by: 

5  cos  a  —  \  n2  —  sin2  a 
-  .  , 

5  cos  a  +  \/q2  ~  sin"  a 

where  5  =  pi/p  (5  =  2)  is  the  ratio  of  densities  and  r/  =  cjc\  (here  q  =  1524/1600)  is  the 
index  of  refraction.  The  reflection  coefficient  at  the  free  surface  is  V2  =  —1. 

Starting  with  upward  and  downward  traveling  plane  waves  emanating  from  a  source  point 
at  an  angle  of  a  from  vertical,  a  series  of  image  plane  waves  are  constructed,  which  when 
combined  with  the  incident  plane  waves,  produces  an  infinite  series  of  plane  waves  satisfying 
the  free  surface  and  interface  boundary  conditions.  Such  a  series  for  the  region  above  the 


T(  a)  = 


2  cos  a 


5  cos  a  +  \Jq 2  —  sin2  a 


(12) 
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source  point  —  Zq  (0  >  z  >  —Zq),  is  found  to  be: 


eHkxx+kvy)  ^2yjyj eikz(z+z0+j2H)  ^_^yn+ny^eikz(z-z0+(j+l)2H) 

3=0 

+  y2eikR~z+zo+32H )  +  (— 1  yn+nyly2eik^~z~zoH3+1)2H) 


(13) 


where  kx  =  sin  a  cos  /?,  ky  =  sin  a  sin  /3,  kz  =  cos  ck.  (Similar  series  for  the  regions  —H  <  z  < 
— zo  and  z  <  —H  can  be  found.) 

Formally  summing  these  series  and  performing  the  integration  in  (3  yields  the  resultant 
integral  for  the  acoustic  held  produced  by  the  (n,  m)  harmonic: 


•m—n  im(f) 
i  C- 


(‘n/2—ioo 


V (a)P™( cos  a)  Jm(kR  sin  a)  sin  ada 


(14) 


where  R  =  \/x2  +  y2  is  the  horizontal  distance  from  the  source  point,  and 


n«)  = 


(i ebz  -  e-bz)(Sz°-H)  +  (_l)n+rab(a)e-J(zo-ff)) 
e~bh  +  Vi(a)ebH 

(— l)"+m(e_bz°  -  ebz°)(Vi(a)eb(z+H}  +  e~b^z+H)) 
e~bh  _|_  V\(a)ebH 


0  >  Z  >  —Zq 


— Zq  >  z  >  —H 


(_1)n+m(e— : bz0  _  ebzo}T(a}e-bi(z+H) 

e~bh  _|_  V\(a)ebH 


—H  >  z 


(15) 


with  b  =  ik  cos  a  and  b\  =  ikyjrj 2  —  sin2  a. 

Evaluation  of  these  integrals  with  slight  modifications  is  detailed  in  Ewing,  Jardetzky, 
and  Press  [9].  For  the  integral  valid  for  0  >  z  >  —  z0,  a  change  of  variables  produces  the 
integral 


where 


2  i 


_  _ _^m—n^im(t> 

k 


Io  =  —I . e 


N(u)  ,  .  nmib 

D(u)Jm^R  ’Pn 


u  du 


Sill  2!/ 


(16) 


N(u) 


<5z/cos[(^o  —  H)v\  +  iv i  sin[(z0  —  H)v]  m  +  n  even 
iSv  sin[(z0  —  H)v\  +  vx  cos[(^0  —  H)v\  m  +  n  odd 


(17) 


D(u)  =  Sv  cos[Hv]  —  iv i  sin  [Hu]  (18) 

v  =  \Jk2  —  u2  ,  Vi  =  \jk\  —  u2  (19) 


This  integral  is  evaluated  using  contour  integration,  by  introducing  a  small  imaginary 
perturbation  of  the  wavenumbers  which  shift  the  poles  of  the  integrand  into  the  first  quadrant 
of  the  complex  u  plane.  A  branch  cut  following  a  hyperbolic  path  emanates  from  the  branch 


5 


point  at  u  —  k\  and  asymptotically  approaches  the  positive  imaginary  axis  (the  integrand  is 
“even”  in  v  so  there  is  no  cut  emanating  from  u  —  k,  and  the  Riemann  surface  is  two-leaved). 
The  Riemann  surface  chosen  is  such  that  the  imaginary  part  of  17  is  positive  (to  ensure  decay 
of  any  waves  propagating  into  the  bottom  half  space). 

The  integral  I 0  is  split  into  two  parts  (integrals  I\  and  J2)  by  rewriting  the  Bessel  function 
Jm(Ru )  =  \[H$  +  H$]  (superscripts  of  the  Hankel  functions  correspond  to  subscripts  of 
the  Fa).  The  path  of  integration  chosen  to  determine  each  of  these  integrals  is  shown  in 
fig.  2.  The  contour  integration  in  the  first  quadrant  correponds  to  the  integral  I±,  while  J2’  s 
contour  is  in  the  fourth  quadrant. 

The  line  integrals  C\  and  Cg  vanish  as  R  — >  00.  The  integrals  C4  +  Cg  and  C5  +  C'10 
cancel  in  pairs  for  the  case  m  is  even  clue  to  the  identity  H$(—Ru)  =  (— 1  )m+1H$(Ru).  In 
the  event  m  is  odd,  the  integrals  can  be  shown  to  decay  exponentially  as  R  — »  00. 

The  branch  cut  integrals  (C'2  and  C'3)  are  combined;  a  change  of  variables  is  performed 
(. x 2  =  kl~u2,  where  along  C'2,  =  —x  and  along  C'3,  17  =  x);  and  the  asymptotic  expansion 

of  for  large  argument  (R)  is  used.  The  combined  branch  cut  integrals  can  be  written: 


with 


F(x,  </>,  z)eiRV^ 


dx 


as  R  — »  00 


(20) 


F(x,<j),z )  = 


2 1 


x 


k  \  ixRi\Jk\  —  x2 
sin[^\/72  —  x2 


eim4,N(Jkl  -  x2)Pnm( 


x z 


x2  sin  [H^/'y2  —  x2 


\/72  —  x2  52( 72  +  x2)  cos2[P\/72  —  x2\  +  x2  sm2[Hy/^2  —  x2 


(21) 


and  7  =  \Jk2  —  k2.  Expanding  F{x,  (f),  z )  in  a  Taylor  series  about  x  —  0,  using  integration  by 
parts,  and  a  stationary  phase  analysis  results  in  a  leading-order  behavior  for  the  branch  cut 
integrals  of  0(1/R2).  Therefore  in  the  far-held  analysis  the  branch  cut  integral  is  negligible 
compared  to  the  contribution  of  the  contour  integral  from  the  poles. 

For  the  present  problem,  there  are  two  real  roots  (up)  to  D(u)  which  lie  in  the  interval 
ki  <  up  <  k.  The  two  roots  are  graphically  displayed  in  fig.  3,  and  have  the  values: 
u  1  ~  1.98522  and  «2  ~  1.93544.  Ultimately  we  obtain  the  far- held  approximation  for  the 
integral  I0: 


h 


27T  i 

Wk 


■m—n  im(f) 
i  C- 


Hup 


P?{vp/k)Hg\Rup)vpH  sin(zvp)S(up) _ 

cos (Hup)  sin (Hup)  —  6~2  sin 2{Hvp)  tan (Hup) 


(22) 


with  equal  to  the  depth  of  the  spherical  shell’s  center,  R  the  horizontal  distance  from  the 
sphere’s  center,  up  =  Jk2  —  w2,  and 


S(up) 


i  sin(^o^p)  m  +  n  even 
cos(^oz/p)  m  +  n  odd 


(23) 


The  residual  calculations  for  the  region  —  z0  >  z  >  —H  yield  identical  results,  whereas 
for  z  <  —H  a  similar  calculation  produces  the  result: 


6 


(24) 


ir-Ueim*  y  P?WQx£,(R'b 

Hkd  p  HiSp  —  cos (Hup)  sin (Hup)  —  5~2  sin 2(Hup)  ta n(Hup) 

Summing  now  over  each  spherical  harmonic,  and  over  each  spherical  shell,  using  the  large 
argument  approximation  for  the  Hankel  function,  and  reordering  the  summations  gives 

N 

E  E  A-mp  (0o)'hmp(0O,^)  (25) 

m=-N  P 


■A-mp^o)  — 


2ni 

Ilk 


vpH 


upH 


e~iUp( Xl0  cos  <t>o+ylo  sin  0o) 
l 


cos (ispH)  sm(upH)  —  5~ 2  snr(ypH)  ta n(vpH) 

l  |  N 

V  Bl  jm~n  pr 

/  v  nm 1  1  n 

n=\m\ 


i  sin  upz0 

COS  VpZft 


m  +  n  even 
m  +  n  odd 


(26) 


'^mp((t>  o) 


gi iRoUp-in/4  ^—m^imtpo 


sin  (upz) 

-sin  (Hup)eiz+H)y/^-^ 


0  >  z  >  -H 
z  <  —H 


(27) 


with  (xjj,  yl0,  Zq)  the  center  of  the  Ith  spherical  shell,  0o  is  the  observation  angle  measured  from 
a  reference  origin  of  the  array,  and  Rq  is  a  reference  far-held  distance  from  the  array  origin. 
Note  that  for  a  general  array,  the  amplitude  has  angular  dependence,  but  for  a  vertical  array, 
values  of  xl0  and  yl0  are  zero,  and  therefore  amplitudes  are  independent  of  the  observation 
angle. 

In  the  problem  studied  here,  a  vertical  array  has  been  used.  The  “eigenfunctions”  for 
the  half  space  problem  (in  the  z  coordinate)  are  shown  in  fig.  4.  Due  to  the  orthogonality 
of  these  eigenfunctions,  when  determining  the  power  output  from  the  array,  the  total  power 
output  is  simply  the  sum  of  the  power  output  from  each  of  the  eigenmodes.  The  total  power 
is  calculated  by  performing  the  integral: 


Rn  {  1  r~H  r 2,r  dv 

lim  — <—  /  /  \p(R0,(/>,z)\\-^-(R0,<f>,z)\  d(j)dz 

Ro^oo  2uj  {  Pi  J-oo  J 0  oRo 

1  r°  r2'K  dp  ) 

+  -/  /  p(-Ro,<M)  k7r(-Ro,<M) \d(j)dz\  (28) 

P  J-H  Jo  oRq  J 

The  power  output  per  mode  for  the  case  of  the  vertical  array  is  given  by: 

n  =  EMrA  | jjU  _  sin (vpH)  cos (vpH)  —  5~2  sin 2(upH)  tan( vpH)\  (29) 

pckupH  L  J 

and  the  total  power  is  then: 


jjtoiai  =  ump 
m,p 


(30) 
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RESULTS 


As  a  first  test  of  the  resulting  formulas,  a  vertical  array  of  two  spherical  shells  are  placed 
in  the  middle  of  the  sound  channel  separated  by  scaled  distances  of  6.0  and  12.4  units.  For 
the  6.0  separation,  the  proportion  of  power  in  the  first  and  second  modes  are  .808  and  .192 
respectively,  while  for  the  12.4  unit  separation,  the  respective  powers  are  1.0  and  0.  Figs. 
5a  and  5b  plot  contour  lines  of  the  downheld  pressure  held  as  a  function  of  depth  and  radial 
distance  in  meters.  As  can  be  seen  from  fig.  5a  (the  6.0  unit  separation),  a  sinusoidally 
shaped  bell  (which  corresponds  to  high  amplitude  or  illuminated  zones)  is  present,  with  a 
wavelength  close  to  that  determined  analytically  from  the  formula  [2]: 

27 r 

A12  = - «  252m  (31) 

u  i  -  u2 

Fig.  5b  with  the  12.4  unit  separation  is  a  case  in  which  power  is  being  transmitted  in 
only  the  first  mode.  There  is  no  apparent  sinusoidal  behavior  in  the  held,  but  one  can  see 
the  effect  of  radial  drop  off  by  the  downward  sloping  pressure  amplitude  contour  lines. 

In  the  next  three  figures,  power  transmission  from  the  two  element  array  are  shown  for 
a  range  of  spacing  between  the  spheres.  In  these  figures  the  shells  are  centered  about  the 
middle  of  the  sound  channel.  Figs.  6  and  7  graph  the  proportion  of  power  transmitted 
into  the  first  and  second  “modes”  of  the  waveguide.  As  can  be  seen,  at  a  separation  of 
about  6 m  all  the  power  is  transmitted  in  the  first  mode  of  propagation,  while  for  decreasing 
separations,  amplitude  of  the  second  mode  increases  to  its  maximum  value  of  about  28%. 
Also  shown  in  these  two  figures  is  the  effect  of  increasing  the  number  of  spherical  harmonics 
included  in  the  code.  The  proportional  power  calculations  are  unaffected  by  an  increase  in 
the  number  of  spherical  harmonics  kept  (. N )  from  2  to  5. 

In  addition,  figs.  6  and  7  highlight  a  curious  feature  of  the  results.  A  third  set  of  displayed 
data  shows  that  the  coupling  between  the  spherical  shells  can  be  discounted  in  determining 
relative  power  outputs.  In  this  third  data  set,  the  radiation  held  amplitudes  ( Blmn )  are  just 
the  values  that  a  single  spherical  shell  would  have  in  an  infinite  fluid  medium.  As  one  can 
see,  the  proportion  of  power  found  by  neglecting  coupling  between  spherical  shells  yields 
nearly  identical  results  to  those  of  the  fully  coupled  case.  Only  a  slight  variation  results 
when  the  shells  are  in  close  proximity. 

If  one  considers  the  total  power  output  from  the  array,  however,  a  distinct  difference 
between  the  independently  radiating  spherical  shells,  and  the  coupled  shells  can  be  seen 
in  fig.  8.  From  this  figure  it  is  apparent  that  the  coupling  between  spheres  can  have  a 
profound  effect  on  the  amount  of  total  power  delivered  by  the  array  -  overestimating  the 
total  power  by  as  much  as  50%  or  underestimating  it  by  as  much  as  20%  when  the  shells  are 
close.  It  should  be  noted  that  the  scaled  wavelength  of  the  radiated  held  is  2n.  Therefore 
the  general  conclusion  drawn  for  this  simple  two  element  array,  is  that  separations  between 
array  elements  which  are  whole  number  multiples  of  the  wavelength  increase ,  while  array 
separations  at  alternate  half  wavelength  decrease  the  total  power  output,  but  do  not  affect 
the  proportion  of  power  distributed  between  propagating  modes  of  the  sound  channel. 

As  a  final  graph,  the  effect  of  adjusting  the  relative  array  origin  as  well  as  the  relative 
separation  of  the  spherical  shells  is  examined.  Fig.  9  is  a  contour  graph  of  the  total  power 
output  from  the  two  element  array  as  a  function  of  scaled  separation  between  array  elements 


and  depth  of  the  array.  (A  line  drawn  horizontally  across  this  graph  is  essentially  the  graph 
in  fig.  8.)  Two  “peaks”  are  apparent.  In  each  of  these,  the  scaled  separation  is  on  the  order 
of  a  wavelength  («  6).  The  peaks  are  also  close  to  either  the  free  surface  or  the  bottom 
interface  of  the  fluid  layer.  This  is  assumed  to  be  due  to  the  fact  that  image  sources  are 
closer  to  the  actual  sources  in  each  case. 


CONCLUSIONS 

This  technique  has  the  advantage  over  strictly  numerical  techniques  in  that  the  problem 
has  been  decoupled  into  manageable  parts.  The  methodology  could  easily  be  extended  to 
consider  actual  transducers  by  the  addition  of  T  matrix  calculations  for  the  given  transducers. 
Such  work  is  in  progress.  Using  finite  elements  to  determine  the  T  matrix  for  a  given  class 
of  transducer,  and  inputting  the  matrix  into  the  above  semi-analytic  technique  can  be  used 
to  determine  an  arbitrary  array  configuration’s  power  output  at  a  prescribed  frequency  of 
operation. 

A  limitation  of  the  above  analysis  is  that  a  relatively  simple  model  of  a  shallow  water 
environment  has  been  employed.  Extensions  to  more  realistic  conditions  might  include  using 
the  radiation  coefficients  found  from  the  coupled  elements  immersed  in  an  infinite  fluid  (Blnm 
from  above)  as  inputs  to  a  parabolic  approximation  method,  or  a  ray  method. 
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FIGURE  CAPTIONS 


Figure  1.  The  Pekeris  model  of  a  shallow  water  channel  of  uniform  depth  containing  an 
array  of  radiating  spherical  steel  shells. 

Figure  2.  The  contour  integral  evaluated  in  determining  the  far  field  response  of  the  radi¬ 
ating  spherical  harmonic  h$(kr)eim<t>P™ (cos  9).  Along  C2  the  real  and  imaginary  parts  of 
v\  are  negative  and  zero,  while  along  C3  z/Rs  real  and  imaginary  parts  are  positive  and  zero. 

Figure  3.  Poles  of  the  Integrand.  Intersections  of  the  curves  tan(V Hu)  and  5v/ivi  as 
functions  of  v. 

Figure  4.  Pressure  Field  Eigenfunctions.  Eigenfunction  profiles  as  functions  of  depth  for 
the  two  trapped  modes. 

Figure  5.  Far-Field  Pressure  Contours  for  two  different  spherical  shell  spacings.  For  the 
top  figure  the  array  center  ias  located  at  the  center  of  the  channel,  and  the  spacing  between 
shells  is  3  meters.  For  the  bottom  figure,  the  array  center  is  again  at  the  channel  center  but 
the  separation  is  now  6.2  meters. 

Figure  6.  Relative  Power  Output  In  First  Propagating  Mode  (n01/IItotai). 

Figure  7.  Relative  Power  Output  In  First  Propagating  Mode  (n02/IItotai). 

Figure  8.  Total  Scaled  Power  Output  (Scaled  by  the  Power  Output  from  a  Single  Spherical 
Shell  in  an  Unbounded  Fluid) 

Figure  9.  Total  Scaled  Power  as  a  Function  of  Shell  Separation  and  Array  Center  Depth 
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